Propagation of thermal excitations in a cluster of vortices in superfluid ^He-B 
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We describe the first measurement on Andreev scattering of tliermai excitations from a vortex 
configuration witli known density, spatial extent, and orientations in "^He-B superfluid. Tlie lieat 
flow from a blackbody radiator in equilibrium rotation at constant angular velocity is measured 
with two quartz tuning fork oscillators. One oscillator creates a controllable density of excitations 
at 0.2 Tc base temperature and the other records the thermal response. The results are compared 
to numerical calculations of ballistic propagation of thermal quasiparticles through a cluster of 
rectilinear vortices. We find good agreement which supports the current understanding of Andreev 
reflection. 

PACS numbers: 67.30hb, 67.57.hc, 67.30.em 
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I. INTRODUCTION 



Measurements of quantized vortices and attempts to 
develop methods for the visualization of different vortex 
configurations has been central in superfluid studies. Re- 
cent interest has focused on quantum turbulence, espe- 
cially in the zero-temperature limitii. In Bose-Einstein- 
condensed atom clouds motion of vortices has been re- 
cently investigated with optical means^. In the long- 
studied case of superfluid "'He, vortices can be detected, 
e.g., by ion trapping on vortex cores or second sound 
attenuation^, while recent work on turbulence studies 
has been making use of transmission measurements of 
charged vortex rings^ or trapping of micron-sized tracer 
particles in vortex tangles^, as well as by analyzing the 
drag force exerted on vibrating structures^. In superfluid 
•^He the traditional method to study vortices is nuclear 
magnetic resonance^. The superfluid flow due to quan- 
tized vortices modifies the order parameter field and thus, 
the NMR signal. In uniform rotation a resolution of a 
single vortex can be obtained^i^ in a measurement of the 
counterflow velocity at temperatures T > 0.5 Tc. At 
very low temperatures, in the limit T/Tc 1, the most 
powerful tool is the Andreev scattering of thermal excita- 
tions. This technique has been developed and exploited 
at the University of Lancaster—. 



Hitherto the Andreev scattering technique has only 
been used to detect turbulent vortex tangles, which for 
interpretation have been assumed to be homogeneous and 
isotropic, but which in practice are of unknown density 
and poorly known spatial extent. Thus, it has not been 
possible to compare theoretical predictions of heat trans- 
port in vortex systems directly to experimental results. 
In this work we provide such a comparison and justify the 
use of the Andreev reflection technique as a visualization 
method of vortices in superfluid ^He-B in the limit of 
vanishing normal fluid density. 



II. ANDREEV REFLECTION FROM VORTEX 
LINES 

In the ballistic regime of quasiparticle transport the 
mean free path of thermal excitations is longer than the 
dimensions of the container. Therefore, thermal equi- 
librium is obtained via interaction of quasiparticles and 
container walls and the collisions between excitations can 
be neglected. In the presence of vortices, the superfluid 
flow field around the vortex lines can constrain the quasi- 
particle trajectories. 

In the rest frame of the superfluid condensate the BCS 
dispersion relation E{p) is symmetrical and the mini- 
mum energy is the pressure dependent superfluid energy 
gap A. The standard picture of Andreev reflection con- 
siders an excitation moving towards an increasing energy 
gai>ii. In •^He-B the superfluid flow fleld modulates the 
minimum in the excitation spectrum. Using the notation 
by Barenghi et alMi, the energy E of the excitation with 
momentum p in the flow field around a vortex is given 
by 
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where ep — /2m* — e^? is the effective kinetic energy 
of the excitation measured with respect to the Fermi 
energy tp and p = |p|. Excitations with > are 
called quasiparticles and excitations with ep < are 
called quasiholes. For quasiparticles the group velocity 
Vg(i?) = dE/dp is parallel to the momentum p whereas 
for quasiholes it is antiparallel. Our experiments are 
performed at the 29 bar pressure, at which the effec- 
tive mass m* « 5.42 7713, where is the mass of a 
■^He atom. The superfluid velocity Vg is proportional to 
the gradient of the phase tp of the order parameter, i.e., 
Vg = fi,/(2?Ti3)V(/3. If we consider a vortex oriented along 
the z axis in cylindrical coordinates {r,(j),z) this becomes 
Vg = k/ (27rr)0, where k = h/2m^ ~ 6.62 x 10~* cm^/s is 
the circulation quantum and (j) the azimuthal unit vector. 

The consequence of the interaction term p • Vg is that 
an excitation traveling at constant energy may not flnd a 
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forward-propagating state due to the superflow gradient 
Vvg along the flight path. When the excitation reaches 
the minimum of the spectrum the group velocity changes 
sign and it retraces its trajectory as an excitation on the 
other side of the minimum. In other words, a quasipar- 
ticle Andreev reflects as a quasihole and vice versa with 
a very small momentum transfer—. 

Let us consider a beam of excitations incident on a sin- 
gle straight vortex. On one side of the vortex the flow 
parallel to the group velocity of the excitation reflects 
quasiparticlcs and on the other side the antiparallcl flow 
reflects quasiholcs. An excitation is Andreev reflected 
if its energy satisfies E < A + pk/ {2Trb) sin 9. Here 9 
is the inclination angle of the excitation trajectory mea- 
sured with respect to the vortex line and b the impact 
parameter. At temperature T the mean excitation en- 
ergy is i? = A + fcsT. In our experiments fc^T ~ O.I A 
and the momentum p is close to the Fermi momentum 
Pf ~ 9.26 • 10^^^ kgm/s, so a typical excitation is re- 
flected if 6 < 5pFKsm6 /{ttA). For an excitation with 
« 7r/4 this translates to ~ Ifim, which is about two 
orders of magnitude larger than the coherence length 
^ sa 15 nm and the vortex core radius. Thus, in a typical 
experimental situation the probability of scattering off a 
vortex core is negligible compared to the cross section for 
Andreev scattering from the flow flcld around the vortex. 



III. EXPERIMENTAL METHODS 

In our experiment we study the heat transported by 
excitations through a cluster of vortices. Bradley and 
coworkers did a similar measurement with a vortex tan- 
gle as the structure reflecting excitations^^. Our experi- 
ment is made in a fused quartz cylinder flllcd with "^He- 
B. The cylinder is divided in two parts with a 0.7 mm 
thick quartz division plate. The lower part consists of a 
30 mm long, 3.5 mm inner diameter tube, which opens to 
a heat exchanger made out of sintered silver. The sinter 
provides good thermal contact with the nuclear cooling 
stage so that the supcrfluid ^He below the 0.3 mm orifice 
can be cooled down to below 0.15 Tq. The upper part 
can be modeled as a blackbody radiator (BBR), an en- 
closure with a weak thermal link to the outside superfluid 
•^He via a small orifice in the division platei^. Our BBR 
consists of a 12 cm long section of the quartz tube with 
6 mm inner diameter. This volume is furnished with two 
mechanical resonators, one acting as a thermometer and 
the other as a heater. The heater is used for generating 
a beam of ballistic quasiparticlcs through the orifice. 

Our resonators arc commercial quartz tuning forks, 
which have recently been characterized for probing quan- 
tum fiuidsi&iii. The forks are made of piezoelectric 
quartz crystals with electrodes deposited on the surface. 
When driven with alternating voltage, the two prongs 
of the fork oscillate in anti-phase producing a current 
/, which is proportional to the prong velocity Vp. The 
heater fork signal is amplified with a room-temperature 
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FIG. 1. The experimental setup. The upper experimental 
volume modeled as a black body radiator is separated from 
the heat exchanger volume at the bottom by a division plate 
with a conical orifice with 0.3 mm diameter. The upper divi- 
sion plate with a 0.75 mm diameter aperture is not relevant 
to the measurement described here. The BBR houses two 
quartz tuning fork oscillators one acting as heater, the other 
as thermometer. 



I/V converter— before being fed to a two phase lock-in 
amplifier. This was found to be important to reduce ca- 
pacitive losses in the signal line, and thus to measure 
accurately the signal amplitude which is proportional to 
the power generated by the fork. The thermometer fork 
has 32 kHz resonance frequency, a prong cross section of 
O.IO mm X 0.24 mm and a length 2.4 mm. The heater 
fork has a higher resonance frequency, 40 kHz, to pre- 
vent any interference between the forks. The prongs 
of the heater are 2.9 mm long and the cross section is 
0.36 mmxO.44 mm. 

In our temperature range the resonance width of the 
tuning fork depends only on the damping from ballistic 
quasiparticlcs. The dependence of the linewidth A/ on 
temperature and prong velocity is given by 



A/ = A/i, 
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where A is a geometrical factor close to unityi^. The sec- 
ond term in the parenthesis is due to Andreev reflection 
of thermal quasiparticles from the potential flow fleld cre- 
ated by the fork prongs moving the liquid around them. 
In our experiments Vp is small and the velocity-dependent 
term in Eq. ([2]) can be neglected. Thus, calibrating 
the fork to act as a thermometer requires determining 
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only the geometry-dependent factor a. The thermome- 
ter is cahbrated at 0.33 Tc against a ■^He-mehing curve 
thermometer, which is thermally coupled to the heat ex- 
changer. Our calibration gives a «17500 Hz for the de- 
tector. The intrinsic damping of the fork was measured 
to be Afint ~ 14 mHz at T ^ 10 mK in vacuum, which 
translates into a quality factor Q ^ 2 • 10^. 

The rough surface of the sinter with a grain size close 
to the vortex core diameter provides excellent spots for 
vortices to nucleate. Thus, the critical rotation velocity 
r^c for vortex formation is lower than 0.1 rad/s in the 
bottom section of the long quartz tube. In the equilib- 
rium vortex state in uniform rotation the sample becomes 
filled with rectilinear vortices oriented along the rotation 
axis. To create the vortex array, which Andreev reflects a 
part of the heat back to the BBR, we rotate our system at 
constant velocity around the axis of the container tube. 
The vortex density in the equilibrium state is determined 
by minimization of the free energy in the rotating frame 
and is given by the solid-body-rotation value n„ = 2^1/ n. 
The array is isolated from the container wall by a narrow 
annular vortex-free layer. The width of the vortex-free 



region ^ K/{\/iVt) is only slightly larger than the inter- 
vortex distance's. 





FIG. 2. Sketch of the experiment. In the system at rest 
(r2=0) all the excitations which do not migrate back to the 
blackbody radiator due to diffuse scattering from the walls are 
thermalized in the heat exchanger at the bottom. In rotation 
(n 7^ 0), part of the beam is Andreev reflected from the cluster 
of vortices. 



To make sure that we have the equilibrium number of 
vortices in the container we first rapidly increase the ro- 
tation velocity to some value which is higher than the 
target velocity for the measurement. Then we go to the 
final velocity and wait for the system to settle to the 
equilibrium vortex state after the annihilation of the ex- 
tra vortices and the slow relaxation of the vortex array. 
In our experimental conditions the relaxation takes about 
one hour. 

All the power entering our experimental volume mod- 



eled as a blackbody radiator must leave through the hole 
at the bottom as a flux of energy-carrying excitations. 
Assuming thermal equilibrium inside the BBR the power 
is given by 

Q{n) ^ j N{E)vg{E)Ef{E)rdEdxdyd(f>d0, (3) 

where N{E) and f{E) are the quasiparticle density of 
states and the Fermi distribution function, respectively. 
In the limit fc^T ^ A the latter reduces to the Boltz- 
mann distribution f{E) = Q-E/ksT ^ rpj^g transmission 
function T = T{E, x, y, 0, 9, fi) is equal to one if an 
excitation leaving the BBR (at position (x, y) on top 
of the orifice to direction ((j>, 0)) reaches the sinter and 
zero if it is reflected back. The integration goes over the 
cross section of the orifice, (j) G (0, 27r), 9 G (0,7r/2) and 
E E (A,oo). The power generated inside the radiator 
can now be expressed as the sum of the Jl-dependent 
residual heat leak Qm to the BBR and the direct power 
Pgen from the excitations produced by the heater fork 
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Te-^iA + kBT)Ah{n). (4) 



Here Afi(^) is the effective area of the orifice, which gets 
smaller when part of the excitations is scattered back 
to the BBR. We omit the flow of excitations from the 
thermal excitation bath in the heat exchanger volume 
since the quasiparticle density there is at least three or- 
ders of magnitude lower than inside the BBR, which is 
at 0.20 Tc. 



IV. STEADY STATE MEASUREMENTS 

In the measurement the heater fork is driven to create 
the desired excitation beam corresponding to power Pgen 
leaving the radiator. By controlling the rotation velocity, 
and thus the vortex density, we can control the fraction 
of Andreev reflected excitations. As illustrated in Fig. [21 
the flow fleld created by the vortices reflects part of the 
beam back to the radiator by Andreev scattering. As a 
consequence, the temperature above the orifice increases 
more than with the same applied heating in the absence 
of vortices. The fraction v of heat reflected back into the 
radiator, which we call the reflection coefficient, can be 
obtained from Eq. (j?]) as 
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There is also a vortex cluster inside the BBR, which may 
cause a small temperature gradient along the cylinder. 
Nevertheless, the main thermal resistance is across the 
orifice and Eq. ([5]) is valid as long as there is thermal 
equilibrium inside the radiator. 

At each rotation velocity, we apply different power in- 
puts to the radiator and measure the corresponding equi- 
librium temperature with the thermometer fork. By plot- 
ting all the temperature-dependent parts in Eq. 1^ as 
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modulating the rotation velocity at different frequencies 
and amplitudes we can study whether the presence of 
these perturbations, which are called Kelvin waves, af- 
fects the reflection coefficient. We find that even an or- 
der of magnitude larger modulation amplitude compared 
to the highest noise peaks in the rotation velocity barely 
affects the fraction of transmitted heat flux. Thus, we 
believe it is safe to omit the effect of Kelvin waves in our 
analysis. Our preliminary measurements at very large 
modulation amplitudes, however, show a decreasing frac- 
tion of transmitted heat, as expected if Kelvin waves are 
generated. In the future, we are hoping to utilize these 
techniques to study Kelvin waves in more detail. 



TRANSIENT MEASUREMENTS 



FIG. 3. Temperature dependent part of the power leaving the 
blackbody radiator as a function of heating power at three dif- 
ferent rotation velocities. The temperature is obtained from 
the linewidth of the detector fork. The shown data points are 
averaged from data measured for about 10 minutes at each 
power. The intercept of the linear fit with the power axis 
gives the residual heat leak to the sample and the effective 
area is given by the inverse of the slope. The slope, the heat 
leak and the scatter in the data all increase with increasing 
angular velocity. The inset shows an example of a detector 
response to a heating pulse starting at time t = 0. 



a function of power Pgen we get a straight line (Fig. . 
From the inverse slope of the line we get the effective area 
Afi and from the intercept with the power axis the heat 
leak Qhi- 

The measurement with no vortices gives Ah{0) ~ 
0.020 mm^. This is less than a half of the geometrical 
area Ag of the orifice. The main reason is the diffusive 
backscattering of the excitations from the walls of the 
0.7 mm thick division plate and the quartz tube below 
it. In any case, the absolute value of the effective area 
is not an important issue since we are only interested 
in the relative change in Eq. ([5|). The heat leak Qm 
varies from 12 pW at ^2 = to 18 pW at fl =1.8 rad/s. 
At high rotation velocities the rotation-induced heat leak 
fluctuates with variations of about 1 pW. The rotation 
velocities used in the measurements had to be carefully 
selected, since mechanical resonances at certain velocities 
cause increased and temporally varying heating. Figure 
m shows the reflection coefficient as a function of the rota- 
tion velocity. In the measured rotation velocity range the 
dependence of ly on the vortex density is approximately 
linear. We believe that the main source of scatter in the 
experimental data comes from the variation of the power 
calibration of the heater fork. 

In the measurements the rotation velocity fluctuates 
on the level of Afi < 0.01 rad/s. Therefore, it is possi- 
ble that we create helical perturbations on vortex lines, 
which can end up increasing the total vortex length and 
decreasing the polarization in our vortex cluster. By 



To test whether our blackbody radiator works as ex- 
pected, we can analyze how the system reaches thermal 
equilibrium when the heater is suddenly switched on. 
The expected time constant for the thermal relaxation 
is T = RC, where the thermal resistance across the ori- 
fice R = (dQ/dT)^^ oc Aj^^ and the heat capacity C is 
given approximately by2i 

C = kBV2^Np (y-^f^ ' (^A + ^^bT^ V. (6) 

Here V « 3.4 cm"^ is the volume of the BBR and Np 
the density of states at the Fermi level. The measured 
time constant is about 25 s (see inset in Fig. [3]), which 
is in a good agreement with the expected time constant 
(~ 32 s) obtained using the effective area from the cal- 
ibration described above. This analysis also proves that 
practically all the heat capacity of the system is in the 
bulk superfiuid '^He. The possible error sources in the 
effective area are the small statistical error in the deter- 
mination of the slope (see Fig. [3]) and uncertainties in the 
power calibration, temperature calibration and the value 
of the gap>22. The reflection coefficient v has a weak log- 
arithmic dependence on the parameter a and depends on 
A only through the temperature calibration [Eq. 
The power calibration, if time- independent, has no effect 
on V. We estimate the overall systematic uncertainty of 
the measurements shown in Fig. 2] to be of the order of 
the scatter in the data. 



VI. SIMULATION CALCULATIONS 

In our numerical simulations we calculate the trans- 
mission function T for our geometry at different rotation 
velocities and solve the integral in Eq. ([3]) numerically us- 
ing Monte Carlo integration with importance sampling. 
The simulations use the exact geometry of our experi- 
mental setup including the thickness and the shape of the 
radiator orifice. Instead of solving for the full equations 
of motions for excitations, which would require too much 
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FIG. 4. The fraction u of the heat Andreev reflected back into 
the blackbody radiator. The temperature inside the radiator 
is 0.20 Tc- The simulation points are obtained by integrating 
Eq. ([S} numerically and solving equations (|4]) and ([S} for v. 



computing power, only the vortices for which the impact 
parameter of the excitation is small enough to allow An- 
dreev reflection are considered. We do not assume perfect 
rctro-rcflection but take into account the small Andreev 
reflection anglei^ Aip = hp~^^ ^ tt / {i£b) . Figure S] shows 
the reflection coefficient as a function of the rotation ve- 
locity from the numerical simulations. The result is in a 
good agreement with the measurements. 

In the simulations diffuse scattering from the container 
walls is assumed based on the results in Ref. [2^ In 
reality, part of the excitations experience Andreev re- 
flection at the container surfaces^!. Including this effect 
would lead to a decrease of the reflection coefficient in 
the numerical simulations, as confirmed by our calcula- 
tions with an unrealistic steplike order-parameter vari- 
ation, and improve the agreement between the simula- 
tions and the measurements. The full accounting for the 
detailed reflection processes at diffusely scattering sur- 
faces, including a slight particle-hole anisotropy^^, would 
complicate the transport calculation and we neglect these 
effects, arguing as follows. Our simulations with fully dif- 
fusive scattering yield a ratio (0) /Ag which agrees well 
with the known Ag and A/i(0) obtained from our mea- 
surement. This is neither the case for the model with 
a steplike order parameter variation at surfaces nor for 
specular scattering, since to produce the measured Ah (0) 



the former would require Ag to be roughly two times 
larger and the latter less than a half of the known geo- 
metrical area. Moreover, quasiclassical calculations with 
a smooth model order-parameter suppression profil o^^i^^ 
show that the real Andreev reflection probability at a dif- 
fusely scattering surface is only a small fraction of that 
in the step model. Finally, even though the details of 
the wall scattering processes do have an effect on 
they do not affect the reflection coefficient v significantly 
compared to other uncertainties in our model. 

Recent numerical studies^liSS indicate that especially 
for dense vortex structures, the total reflecting "Andreev 
shadow" is not necessarily the sum of shadows of single 
vortices. Our clusters are relatively sparse and moreover, 
after the first diffuse scattering from the walls the proba- 
bility for the excitation to migrate back to the radiator is 
not sensitive to small changes in its trajectory. Thus, we 
believe that our somewhat simplified model reproduces 
the real experimental situation with good accuracy. The 
simulations were tested at different hole radii and po- 
sitions of the hole on the division plate. We find that 
the largest effect on the final result comes from the un- 
certainty of the radius: increasing or decreasing it by 
50% changes the reflection coefficient by approximately 
±10%. 



VII. CONCLUSIONS 

In conclusion, we describe the first measurement of 
the interaction between thermal excitations and quan- 
tized vortices in a well-defined configuration. Numeri- 
cal simulations reproduce the experimental results within 
the margin of uncertainty in the model. Quasiparticle 
beam techniques are currently the most popular measur- 
ing method of vortices in ^He-B below 0.2 Tq. Our work 
provides a rigorous quantitative basis for their use and 
further development for direct visualization purposes. 
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